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Abstract. Based on general arguments the in-medium quark propagator in a quark-gluon plasma 
leads to a quark dispersion relation consisting of two branches, of which one exhibits a minimum 
at some finite momentum. This results in a vanishing group velocity for collective quark modes. 
Important quantities such as the production rate of low mass lepton pairs and mesonic correlators 
depend inversely on this group velocity. Therefore these quantities, which follow from self energy 
diagrams containing a quark loop, are strongly affected by Van Hove singularities (peaks and gaps). 
If these sharp structures could be observed in relativistic heavy-ion collisions it would reveal the 
physical picture of the QGP as a gas of quasiparticles. 
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1. Introduction 

Recently, Van Hove singularities in the QGP [1-6], caused by the collective quark modes 
appearing due to the interaction with the medium, have attracted great attention as they 
could probe the QGP as a gas of quasiparticles formed in heavy-ion collisions. One of the 
quark modes, the plasmino branch, has a minimum at some finite momentum leading to 
a vanishing group velocity. Thus the density of states, which is inversely proportional to 
the group velocity, diverges, causing Van Hove singularities on relevant quantities like the 
dilepton production rate and mesonic correlators. Here we would like to discuss the role 
of Van Hove singularities in the QGP. 

In the next section, we briefly recall the classic idea of Van Hove singularities in the den- 
sity of states as discussed in the context of condensed matter physics. We will demonstrate 
in the following two sections that Van Hove singularities appear as a dynamical aspect 
of QGP on account of the effective quark propagator containing the quark self-energy. In 
sec. 3 we will demonstrate this through the hadronic correlators using the hard thermal loop 
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(HTL) resummed qurak propagator whereas in sec. 4 through the nonperturbative dilep- 
ton production rate, using an effective quark propagator obtained by taking into account the 
nonperturbative gluon condensate above the critical temperature derived from lattice QCD. 
In sec. 5 we will argue that the minimum in the plasmino branch and thus the appearance 
of Van Hove singularities is a general property of massless fermions at finite temperature. 
Finally, in sec. 6 we will discuss the possibilities and difficulties to observe Van Hove 
singularities which could reflect the physical picture of the QGP as a gas of quasiparticles 
produced in heavy-ion collisions. 

2. Van Hove Singularities 

In condensed matter physics Van Hove singularities in the density of states were discussed 
on structural aspects of a solid by Van Hove [7] in 1953. The density of state of a system 
{e.g. phonons, electrons) is given as 



implying the number of allowed wave vectors k in the level n having energy a;„(fe) within 
the energy interval u and u + du>. By surface integration it can be given [8] as 



which is explicitly related to the group velocity, | Va;„(fc)|, which is the derivative of the 
dispersion relation. Due to symmetries in a crystal the group velocity vanishes at certain 
momenta, resulting in a divergent integrand in (2). This divergence is integrable in 3- 
dimensions, leading to a finite density of states. In lower dimensions, however. Van Hove 
singularities appear. For example, a 2-dimensional electron gas shows logarithmic singu- 
larities, which have been discussed in connection with high-Tc superconductors [9]. A 
2-dimensional electron gas has also been discussed widely in the context of the quantum 
Hall effect. As discussed here, the Van Hove singularities are linked to the structural aspect 
of a solid whereas, as we will see below, in the QGP they are due to dynamical reasons. 

3. Thermal Hadronic Correlaton Function 

3. 1 Definition 

Meson correlators are constructed from meson currents Jm{t.x) = q{T. x)Tm<i{t,x), 
where Tm = 1, 75, 1^, 7/i75 for scalar, pseudo-scalar, vector and pseudo-vector channels, 
respectively. The thermal two-point functions in coordinate space, Gm{t, x), are defined 
as 




n 



(1) 




(2) 



Gm{t,x) = {JM{T,x)jl{0,d)) 




(3) 
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where t G [0, 1/T], and the Fourier transformed correlation function XM{^n,p) is given 
at the discrete Matsubara modes, a;„ = 2mrT. The imaginary part of the momentum space 
correlator gives the spectral function (jm{<^,p), 

XM{(^n,P) = - / du! . ""^(^'^^ =4. (jm{oj,P) = -Im xm(w,p). (4) 

Using eqs. (3) and (4) we obtain the spectral representation of the thermal correlation 
function in coordinate space at fixed momentum {(3=1 /T), 

Gm{t,p)= du>aM{oJ,p) ■ . . (5) 



3.2 Hard Thermal Loop (HTL) Approximation 

Perturbative QCD at finite temperature is based on the fact that the temperature dependent 
running coupling constant is small at high temperatures due to asymptotic freedom, T — > 
00 =J> as(T) = ,g^/47r ^0. At a typical temperature of T = 250 MeV we expect 
as = 0.3 - 0.5. This suggests that perturbative theory could work at least qualitatively. 

However, restricting only to bare propagators and vertices perturbative QCD can lead to 
serious problems, i.e., infrared divergent and gauge dependent results for physical quanti- 
ties, e.g., in the case of the damping rate of a long wave, collective gluon mode in the QGP. 
The sign and magnitude of the gluon damping rate was found to be strongly gauge depen- 
dent [10, 1 1]. The reason for such behaviour is the fact that bare perturbative QCD at finite 
temperature is incomplete, i.e., higher order diagrams missing in bare perturbation theory 
can contribute to lower order in coupling constant. In order to overcome these problems, 
the HTL resummation technique, an improved perturbation theory, has been suggested by 
Braaten and Pisarski [12] and also by Frenkel and Taylor [13], in which those diagrams 
can be taken into account by resummation. The following prescription has been suggested 
for the improved perturbation theory: 

1. Isolate those diagrams which should be resuimned. The starting point is the separa- 
tion of scales in the weak coupling Hmit, g « 1, since there are two momentum 
scales in a plasma of massless particles: (i) hard, where the momentum ^ tempera- 
ture, T and (ii) soft, where the momentum ~ thermal mass (~ gT, g « 1). These 
diagrams are HTL self energies and vertices given by one loop diagram, where the 
external momenta are soft and the loop momenta are hard. For example, in the case 
of self energies in QED or QCD a result proportional to g^T^ was obtained, which 
is equivalent to the high temperature approximation, T >> p, \po\, found earlier 
by Klimov [14] and Weldon [15], and also to the semiclassical approximation by 
Silin [16]. 

2. Effective propagators in the HTL approximation are constructed using the Dyson- 
Sch winger equations. The HTL vertex follows simply by adding the HTL correction 
to the bare vertex. This effective vertex is also related, except for scalar particles, to 
the HTL propagator through Ward identity. 
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3. These effective quantities, propagators and vertices, can then be used as in ordinary 
perturbation theory if all legs are soft, leading to results that are complete to leading 
order in the coupling constant, gauge independent and improved in the infrared be- 
haviour If, on the other hand, at least one leg of the propagtors or vertices is hard, 
bare quantities are sufficient. 



3.3 Quark self energy in HTL-approximation 



Figure 1. Quark self energy having quark momentum K. 

The most general ansatz for the fermionic self energy in rest frame of the plasma is given 
by [17] 

E{K) = -a{ko,k)If-b{ko,k)jo, (6) 

where K = (fco , k), /c = |k| and the quark mass is neglected assuming that the temperature 
is much larger than the quark mass, which holds at least for u and d quarks. The scalar 
quantities a and b are given by the traces over self energies in figure 1 as, respectively, 

«(^o, k) = ^ [tr (-^S) - fco tr (70S)] , 

6(fco, k)^^ [K^ tr (70S) - fco tr (^E)] , (7) 
and obtained in HTL-approximation as 



a(fco,fc) = — (^l-_ln^j 



biko,k) = ^ ( -fco + In^ ) , (8) 



fcp k^ ] ^0 
fc2 V ' 2fc '^fc7 

where = is the effective quark mass. The quark self energy in eq.(6) has an 

imaginary part below the light cone, fcg — fc^ < 0, representing Landau damping for space 
like quark momenta. Furthermore, the general ansatz in eq.(6) is also chirally invariant in 
spite of the appearance of an effective quark mass. 



3.4 Effective Propagators and Vertices in HTL-approximation 

Resumming the quark self-energy by using the Dyson-Schwinger equation, the effective 
quark propagator in figure 2 can be written as 

s{k)^[^-j:{k)]-\ (9) 
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Figure 2. Effective quark propagator in HTL-approximation 

For massless quarks, decomposing the effective quark propagator in (9) into helicity eigen- 
states, we get 



S{K) = 



2D+{K) 2D_{K) 



with 



Z?±(fco, fc) = (1 + a(fco, k)) {~ko ±k) + b{ko, k) . 
It QED hke HTL-vertex is related to this pro 



(10) 



(11) 



The relevant QED like HTL-vertex is related to this propagator through the Ward iden- 
tity [1,13] 

(12) 




Figure 3. Quark dispersion relation in HTL-approximation along with the free 
(u! — k) one. 

Now, the zeros of D±{K) give the in-medium dispersion relation for quarks. As shown 
in figure 3 the upper curve u!+{k) corresponds to the solution of D+{K) = 0, whereas 
the lower curve uj^ (k) represents the solution of Z)_ {K) — 0. Both branches start from a 
common effective mass, m^, obtained in the k ^ limit [1]. The uj^{k) branch describes 
the propagation of an ordinary quark with thermal mass, and the ratio of its chirality to he- 
licity is +1. On the other hand, the uj^ (k) branch coiTesponds to the propagation of a quark 
mode with a negative chirality to helicity ratio. This branch represents the plasmino mode 
which is absent in the vacuum but appears as consequence of the medium due to the broken 
Lorentz invariance, and has a shallow minimum. This corresponds to a purely collective 
long wave-length mode, whose spectral strength decreases exponentially at high momenta. 
For high momenta, however, both branches approach the free dispersion relation. 
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3.5 Hadronic spectral and correlation function in HTL-approximation 

The hadronic spectral functions, f7M(w,p) in (4), of the temporal correlators are propor- 
tional to the imaginary part of the quark loop diagram. Here we will calculate the imaginary 
part of the quark loop diagram corresponding to the vector meson channel, which is related 
to the static dilepton production rate (p = 0) as [1] 

a.H = -Imx.H = ^ (e^- - ^W^^^i^= 0)- (13) 

For the details of the other channels, e.g., pseudovector, scalar and pseudoscalar, see 
Refs. [5,6]. Now the vector meson correlation function can be written from figure 4 as 







Q = P-K 

Figure 4. Self-energy diagram involving quarks loop in HTL-approximation. 



Xv{uj,p) = 2NcTY, J j^tT[S{K)r^{-Q,-K;P)S{Q)r''iQ,K;-P)] , (14) 

ho 



where, Q = P — K, Nc = 3 and xv denotes the trace over the Lorentz indices of Xm- 
Using the effective quantities from (10) and (12), eq. (13) at p = can be written as 

ay(w) = ^(e'''^ - 1) / dkk^ dx dx'np {x)nF (x) 6 {co - x - x') 

TT Jo J-oo J-oo 



x^+x'^-2e-2ml\' , . . , . 
xp+ [x, k)p+{x ,k) + \1 — ^ I p- {x, k) p-{x , k) 



(15) 



where x and x' are the energies of the internal particles, n p is the fermi distribution func- 
tion, and p± are the spectral function for collective quark modes in the medium [1] having 
a pole contribution from the in-medium dispersion relation and a cut contribution for space 
like quark momenta related to Landau damping below the light cone (fcg < k^). The meson 
spectral function, constructed from two quark propagators, will have pole-pole, pole-cut 
and cut-cut contributions, crM(w) = crPP(a;) -|- c7P'^(a;) -|- a^^{u:i). Here we quote only the 
pole-pole expression for vector mesons. 
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n-PP - 



M.(--(fc2)) 1 2a;Ufc.) 



2kl np (w+ (fcs)) {-uj- (fca)) 



(fcg) _ kl) {col (fcg) - fcf ) 



^ ljI (fca) + jks) - 2fc| - 2ml ^ 



2k3iO 



(16) 



where fci is the solution of w — 2u)^{k) = 0, and ^2 and kl are the solutions of 
UJ — 2ti)_(fc) = and uj — u}+{k) + w_(fc) = 0, respectively, which can have two so- 
lutions at low momenta due to the plasmino mode. Each term in (16) corresponds to 
different physical processes [1,5,6] and is inversely proportional to the derivative of the 
dispersion relation (group velocity). The other contributions (cr^'^ and a^^) can easily be 
obtained [1,5,6] from (15). Similar results for the scalar channel have also been derived 
and detailed discussions are given in Ref. [5,6]. Free meson spectral functions are also 
obtained using bare propagators and vertices in figure 4 as [5,18] 



^m"" (w) = tanh ttM 



(17) 



where qm = (+) — 1 for (pseudo)scalar and (— ) + 2 for (pseudo)vector. 



1E-01 



1E-04 
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cut-cut i-r 

oj/T 



Figure 5. (a) Pseudoscalar and (b) vector meson spectral functions in HTL approxi- 
mation for rUq/T = 1. The free one is represented by crosses. 



In figure 5 we show the different contributions for pseudoscalar and vector mesons for 
niq/T = 1, i.e., g = \/6. Generically the pole-pole term of the meson spectral function in 
HTL-approximation describes three physical processes: (1) the annihilation of collective 
quarks (e.g., P* term in (16)), (2) the annihilation of two plasminos (2"^ term in (16)) 
and (3) the transition from upper to lower branch (3'''^ term in (16)). The transition process 
starts at zero energy and continues until the maximum difference uj — 0.47 rriq between the 
two branches at fca = 1.18 niq. At this point a Van Hove singularity is encountered due to 
the vanishing denominator (^^(fcs) — <^'- (fcs) = 0, implying a diverging density of states. 
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The plasmino annihilation starts at a; = 1.86 m.q with another Van Hove singularity cor- 
responding to the minimum of the plasmino branch at = 0.41 uiq, where ui'_{k2) = 0. 
This contribution falls off rapidly due to the exponentially suppressed spectral strength 
of the plasmino mode for large energies, where only the first process, quark-antiquark 
annihilation starting at w = 2mg, contributes. For large energies this dominates and ap- 
proaches the free results (crosses) for a; >> rriq. The pole-cut and cut-cut contributions, 
which involve external gluons as can be seen by cutting the HTL quark self energy, lead 
to a smooth contribution to the spectral function. The pole-pole and pole-cut contribu- 
tions in both channels are of similar magnitude, while the cut-cut contribution at small lo 
vanishes in the pseudoscalar channel but diverges in the vector channel. This has its ori- 
gin in the structure of HTL quark-meson vertex [1,13], which is required for the vector 
channel containing a collinear singularity [191, whereas bare vertices are sufficient for the 
(pseudo)scalar channel. This implies that higher order diagrams in the HTL expansion 
contribute to the same order in the coupUng constant [20], which, of course, indicates that 
the low frequency part of the vector meson spectral function is inherently nonperturbative. 
As discussed, the vector meson spectral function is linked to the dilepton production [1,5] 
at high temperature whereas the pseudoscalar correlator is related to the chiral condensate 
and thus to the chiral susceptibility [21,22]. 




The temporal correlators for pseudoscalar and vector mesons can be obtained from their 
spectral function according to (5). The competing effects of pole and cut contributions 
in the spectral function carries over to the correlation functions, which is found to be al- 
most similar [5,6], unlike in lattice calculations [23], to the free one with bare propagators 
and vertices. A linear divergence, due to the effective HTL-vertices, of the spectral func- 
tion in the vector channel at low frequencies also renders the temporal correlator infrared 
divergent. Although the pseudoscalar correlation functions are infrared finite, the low fre- 
quency part of the spectral functions will also be modified significantly from contribu- 
tions of higher order diagrams. It is reasonable to consider modified correlation functions, 
A(S'm(t) = Gm{t) — Gm{PI'^), which are less sensitive to details of the low frequency 
part of the spectral functions. In the subtracted correlation functions the infrared diver- 
gences are eliminated. As shown in figure 6 they are well-defined in both channels. One 
observes that after elimination of the infrared divergent parts the structure of the pole and 
cut contributions is similar in the scalar and vector channel. The vector correlator seems to 
be even closer to the leading order perturbative (free) correlator than the scalar one. 
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4. Gluon condensate and nonperturbative dilepton production 



It seems that quite generally nonperturbative effects [24] dominate the temperature regime 
attainable in heavy-ion collisions and the application of perturbative calculations there be- 
comes questionable as the coupling constant is not small. Such nonperturbative effects 
were made expUcit by lattice QCD calculations [23] below and above the phase transition 
temperature. The gluon condensate, which appears due to the broken scale invariance at 
zero as well as at finite temperature, describes the nonperturbative imprints of the QCD 
ground states and have extensively been used for studying hadron properties at zero and 
finite temperature. Lattice QCD so far is not capable to compute dynamical quantities. As 
an alternative method to lattice and perturbative QCD it has been suggested to include the 
gluon condensate obtained from lattice QCD into the parton propagators at zero [25] as 
well at finite [26] temperature. In this way, nonperturbative effects are taken into account 
within the Green functions technique, which then can be exploited to calculate various 
physical quantities Uke the dilepton rate, hadronic correlators, susceptibiUties, etc. 



4.1 Quark propagation 



Figure 7. Quark self-energy diagram in the presence of the gluon condensate. 

Now, the functions a{ko, k) and b{ko, k) in the general effective quark propagator at 
finite temperature defined in (10) and (11) can be obtained using self energy diagram with 
the gluon condensate in figure 7 as [2,3,26] 



9l± 
6 

4 



{£')t -ilk'- kl 



{B')t 



9^ 



kl{£% 



(18) 



with the chromoelectric and and chromomagnetic condensates expressed [2,3,26] in terms 
of the space like (A^) and time like (A,-) plaquette expectation values as 



IT 

as 



11' 



-T^A, 



2 
11 



TT 11 



2 

TT 



(19) 



where (G^)t=o = (2.5 ± 1.0)T^ and the plaqutte values are taken from lattice calcula- 
tions [24]. This propagator is also related to an effective vertex [27] through the QED like 
Ward identity. The dispersion relation in the presence of the gluon condensate is shown 
in figure 8, which shows the same qualitative features as in the HTL case (figure 3) al- 
though it is a completely different approximation. However, the thermal mass of the quark 
is rrieff ~ 1.1 5T, which is independent of the coupling constant unlike in the HTL case. 
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l/Tc l/Tc 
Figure 8. Quark dispersion relation in the presence of the gluon condensate with 
T = l.lTc (left) and T = AT^ (right). Here, read l=fe. 

4.2 Nonperturbative dilepton production 




Figure 9. Dilepton production in the presence of the gluon condensate for T = l.lTc, 

2Tc and 4Tc. 

The dilepton production rate using the gluon condensate can be obtained from the imag- 
inary part of the photon self energy related to (13) 

d^xdvd^p " 6^ JPe^/T_i (20) 

where E = uj = is the energy of the virtual photon with invariant mass M and 

momentum p, and a is the fine structure constant. The nonperturbative dilepton production 
rate has been obtained in Ref. [2,3] neglecting an effective quark-photon vertex containing 
the gluon condensate. The results are shown in figure 9 for p = (left) and p = Tc (right) 
for different temperatures. Asinthecaseof softdileptons[l] and hadronic correlators [5,6] 
calculated within the HTL approximation, we also find peaks and gaps in the dilepton rate. 
The peaks (Van Hove singularities) are caused again by the presence of the minimum 
in the plasmino dispersion (figure 8). The contribution at small M ending with a Van 
Hove singularity comes from an electromagnetic transition from the upper branch to the 
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lower branch ((7+ — > 7*5-) of the quark dispersion relation. After this peak there is a 
gap before the channel for plasmino annihilation q_ 7*) opens up with another 
singularity. This contribution drops quickly but at M = 2mca the contribution from the 
annihilation of collective quarks (<?+ g+ 7*) sets in. This contribution dominates and 
approaches the one coming from the annihilation of bare quarks (Born term [28]) at large 
M. The peaks and gaps are smoothed out to some extent for nonzero p as shown in the right 
panel. It should be noted that there are no smooth cut contributions in contrast to the HTL 
dilepton rate, because the quark self energy containing the gluon condensate (figure 7) has 
no imaginary part. 

5. General quark dispersion relation 

It has been demonstrated in the preceding two sections with two completely different ap- 
proximations that the in-medium dispersion relation reveals identical qualitative features: 
(1) there are always two branches; (2) both of them lie above the free dispersion curve 
(u) = k) and start at the same energy, i.e., effective mass, in the isotropic limit (fc — > 0); (3) 
the two branches have opposite slopes in the isotropic limit; (4) as shown in the two inde- 
pendent examples, both branches approach the free dispersion relation thus the plasmino 
branch must always have a minimum. Considering the most general ansatz for the fermion 
self energy and thus for the quark propagator, it has been argued [4] that the plasmino min- 
imum is a general feature responsible for the appearance of the van Hove singularities in 
the dilepton production rate and hadronic correlators. 

6. Discussions and conclusion 

Considering two independent approximations, the HTL and the gluon condensate, for the 
quark propagator, we have studied the in-medium quark dispersion relation. The disper- 
sion relation, in both approximations, is found to govern general characteristics implying 
the physical picture of the QGP as a gas of quasiparticles. The important feature, which 
appears to be general, is the minimum in plasmino branch due to the broken Lorentz in- 
variance. 

One loop dispersive curves are gauge invariant and whether they are physically observ- 
able is, again, a highly nontrivial question. Future lattice calculation might be able to 
investigate the dispersion relation [29]. However, it is easy to measure explicitly photon 
dispersive curves in a QED plasma by studying the propagation of classical electromag- 
netic waves, which measures the electric charge density as a function of time and spatial 
coordinates and thus determine the frequency and wave vectors. 

But due to confinement there is no such thing in nature like color fields and also device 
which could measure the color charge density. What one can measure are correlators of 
colorless currents Uke mesonic correlation functions and the dilepton production rate. As 
seen, they will be affected by the in-medium dispersion relation because a colorless current 
is always coupled to a pair of colored particles through their Green's functions which are 
related to the characteristics of the dispersion relation only in a indirect way. Also, the 
thermal modification of vertices is important as they are related to the thermal Green's 
functions through the Ward identity. 
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We have calculated the dilepton rates and hadronic correlators using the in-medium 
dispersion relation obtained in two different approximation and have found in both cases 
Van Hove singularities, which are due to the general characteristic (plasmino mode) of 
the dispersion relation. However, to some extent these structures will be smeared out by 
consideration of finite momentum, higher order processes and the space-time evolution 
of the fireball. If these structures are observed in the dilepton spectrum or in hadronic 
correlation functions in relativistic heavy-ion collisions, they would certainly reveal the 
general nature of the dispersion relation implying the physical picture of the QGP as a gas 
of quasiparticles. Such a signal is not expected from the hadronic phase [30]. These effects 
are interesting topics of future investigations. 
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